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Abstract 



We consider a simple model for the growth of isolated steps on a vicinal 
crystal surface. It incorporates diffusion and drift of adatoms on the terrace, 
and strong step and kink edge barriers. Using a combination of analytic meth- 
ods and Monte Carlo simulations, we study the morphology of growing steps 
in detail. In particular, under typical Molecular Beam Epitaxy conditions 
the step morphology is linearly unstable in the model and develops fingers 
separated by deep cracks. The vertical roughness of the step grows linearly in 
time, while horizontally the fingers coarsen proportional to t^'^^. We develop 
scaling arguments to study the saturation of the ledge morphology for a finite 
width and length of the terrace. 

PACS numbers: 68.55.Bd, 68.35.Fx, 82.20.Wt 
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I. INTRODUCTION 



Atomistically controlled growth of metal and semiconductor crystal surfaces constitutes 
an important field of research both from technological [|l| and fundamental theoretical 
points of view. Among all the different growth techniques, Molecular Beam Epitaxy (MBE) 
has a special status since it can be very efficiently used to produce growth in well defined 
layer-by- layer growth mode. Experiments using the reflection high-energy electron diffrac- 
tion technique indicate two main mechanisms of growth in such cases: layer growth by 
nucleation and spreading of 2D islands on a nominally fiat substrate, and step-flow growth 
of a substrate with steps. In the latter case, it is crucial to be in the regime where the flux 
of adatoms is small enough, and their diffusion fast enough to avoid island nucleation on 
terraces. Such a window of the relevant physical parameters may be found experimentally 
for many materials 

An important practical realization of the step-flow situation is epitaxial growth on a 
vicinal surface that is cut in a direction slightly off from a high-symmetry one. Such surfaces 
often consist of broad terraces of size separated by monoatomic steps at distance i from each 
other. By changing the miscut angle, the density of the steps and thus i may be controlled. 
The physics of MBE growth on such surfaces can be in the simplest terms described by 
the following schematic model (see Fig. 1). There is a flux F of adatoms that impinges 
upon the terraces. Particles on terraces then diffuse around with an associated diffusion 
constant D, and may be desorbed after a time r. Upon approaching step edges, particles 
can either cross them from above or below, be reflected from them, or be incorporated into 
the growing ledge. Attachment is usually characterized by Arrhenius type rate constants 
k+ and A;_ which refer to the average rates of particles arriving at the ledge from below or 
above, respectively. 

This simplified picture of step flow growth was first introduced by Burton et al. . More 
recently, attention has been drawn to the fact that in many real systems, k+ and /c_ need not 
be equal because of the existence of step edge harriers . These barriers may often be 
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present at step edges due to reduced coordination of atoms. Recent theoretical work shows 
that the step barriers play an important role in controlling growth under MBE situations 
In particular, if these barriers are high adatoms cannot cross steps, and the particle 
current will be in the direction of ascending steps. For vicinal surfaces, this stabilizes the 
step flow growth mode when nucleation on terraces is neglected. If the average distance 
between nucleation centers is ^at, step flow growth requires that l/dN ^ 1- 

Most of the recent work dealing with step growth has concentrated on the global prop- 
erties and kinetic roughening of growing surfaces with steps under MBE conditions 
However, attention has also been paid on the properties of individual steps under growth 
10|-|T4[. It is a well known property of ideal, isolated ID steps that they are thermally rough 



above zero temperature due to kinks. Using linear stability analysis. Bales and Zangwill |TT 
have shown that in a system with unequal attachment rates ^ k_, a straight terrace 
ledge can become unstable when > This kind of growth-driven instability is par- 
ticularly interesting since it may lead to the appearance of "wavy" patterns of the ledges. 
More recently, Salditt and Spohn have argued that in addition to the instability, there 



is a regime for isolated steps (in the case of strong step edge barriers) where the Kardar- 
Parisi-Zhang (KPZ) theory of kinetic roughening is valid. In this regime, the width of 
the ledge eventually roughens as t^^^ in analogy to many ID surface deposition models f^. 

In this work, we study the nature of ledge or step edge morphologies, and the question of 
their roughening behavior in a simple but nontrivial model of isolated steps. This model is in 
part motivated by the energetics of adatoms on Si(OOl) surfaces with widely separated steps. 
In the model, we assume infinitely strong step edge barriers, and biased diffusion both on the 
terrace and at the ledge. As expected |Tl| , the ledge always becomes morphologically linearly 



unstable due to the dominance of the one-sided diffusion field. Through a combination of 
analytic arguments and computer simulations we show that the ledge develops finger-like 
structures and its roughness grows linearly in time, in contrast to the KPZ type of roughening 
predicted by Salditt and Spohn |]T2[ in the stable regime. In addition, we study the lateral 
coarsening of these fingers and show that it follows a t^'^^ behavior. We develop scaling 



arguments to study the influence of tlie finite width and finite length of the terrace on the 
growth. Finally, we discuss the relevance of these results with respect to experiments on 
ledge roughening under MBE growth [jl6 . 



II. ANISOTROPIC STEP GROWTH MODEL 
A. Definition of the Model 

The model is defined on a two dimensional square lattice where there is a single growing 
step. The average direction of the ledge is along the x axis where the width is with 
periodic boundary conditions. Initially at time t = the step at ?/ = is completely straight 
with no thermal fluctuations present. Growth of the ledge is initiated by depositing a single 
particle randomly on an empty, randomly chosen terrace site in front of the ledge at ?/ > 0. 
After this, the particle performs random walk and drifts towards the ledge by jumping id 
lattice sites in the —y direction at every random walk step on the average. This means that 
during each step, the particle moves in the —y direction with a probability (l/4+£rf)/(l+£rf), 
while for the other three directions the probability is (l/4)/(l + f^). 

The ledge acts as an absorbing boundary to the particle with the following rules (see 
Fig. 2): (i) if the particle arrives at the "top" (a section along the x direction of the step), 
it is incorporated into it and becomes immobile; (ii) if the particle arrives at the "side" (a 
section along the y direction of the step), it will instantaneously slide down along the ledge 
to the —y direction until it reaches the corner site where it is permanently incorporated into 
the step. These rules guarantee that the set of step heights {h{x, t)} as measured from y = 
obey the solid-on-solid restriction, and the step forms a compact structure. 

After the particle has been incorporated into the step, a new particle is deposited and 
the process is repeated. Time in the model is measured in terms of the average height 
of the growing step edge. We note that the size of the terrace in the y direction is not 
fixed, but is chosen in such a way that the distance from the highest point of the step 
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H(t) = max {h{x,t)} is kept at a constant value. The corresponding boundary above 
y = H + Ly + 1 is completely reflecting and remains straight. This means that a particle at 
y = H + Ly that takes a step in the y direction, is immediately reflected back. 

An important feature of the growth model is the deposition of adatoms on the lower 
terrace only. This is tantamount to assuming that the step barriers are inflnitely high with 
A;_ = so that adatoms are reflected from a downward step leading to an average particle 
current in the —y direction towards the up steps. The drift term id is defined only in an 
average macroscopic sense and will depend on the deposition flux and the concentration of 
adatoms on the terrace in front of the step. Also, since we assume that there is no desorption 
of adatoms (r = oo), i^. also depends on the velocity of the step which in turn depends on 
the terrace length Ly. Thus i^. is, in principle, determined self-consistently by the other 
parameters of the model but we regard it as an independent parameter which may be varied 
externally [|17 . 



Finally, we would like to mention that the growth rules of the model are in part motivated 
by adatom dynamics on Si(OOl) surfaces with widely separated steps [P^S|-pT[] under typical 



MBE conditions. Namely, on Si(OOl) diffusion is spatially anisotropic both on the terrace 
TB|-pD| and at the step edges |TB|,|TP||. However, at least for the case of single- height 



steps on Si(OOl), microscopic calculations ||l8i|l9i and experiments |]T6i^ indicate that there 



is no signiflcant step edge barrier. Thus, we make no attempt to realistically model the 
complicated adatom dynamics in this system, since the main motivation here is to study the 
generic features of the unstable regime for isolated steps. 

B. Simulation Algorithm for the Model 

A straightforward Monte Carlo simulation of the growth model proposed here is in prin- 
ciple possible, but very difficult for large values of Ly and small drifts. This is because 
particles landing on the terrace may wander arbitrarily far from the step edge, and thus 
the time for a particle to become incorporated into a growing step may become very large. 
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This problem can be solved by considering the properties of 2D random walkers on a finite 
or semi-infinite plane. For such cases, it is possible to calculate analytically the spatial and 
temporal probability distributions for the walkers. The idea then is that for particles that 
initially land on the terrace with y > H + 1 (which is always the case if Ly — oo), the 
simulation can be started by releasing them from an imaginary line that runs along the x 
direction just one lattice site above the highest step, i.e. at y = H + 1 (see Fig. 2). If 
the particle crosses the line again in the +y direction while performing random walk, it 
is immediately returned to it with a new x coordinate chosen from the appropriate spatial 
distribution which will be derived below. In the Appendix we also calculate the mean arrival 
(first passage) time of a walker and indeed show that this time becomes very large for small 
values of id and large Ly. 

More specifically, to implement the simulation algorithm described above, we need to 
calculate the spatial probability distribution function PLx,Ly{x) which is used to obtain the 
new position for a walker that crosses the line y — H + 1 at any point. In other words, 
a walker crossing the line being at {xo,H + 2) with any xq, is put back to the new site 
{x — Xq, H + 1) with the probability Piyix) where we assume for simplicity that — oo 
(see Fig. 2). For a discrete walker, this function satisfies the recursion relation 




7^[Wx,0 + PLyix - 1) + PLy{x + l) + 



oo 



Yl PLy{x - y)PLy-l{y)], 



(1) 



where a = 2 + 21^ and b 



1 + Md. Using the standard Fourier transformation 



oo 



hyik) = E ^''''PLyix), 



(2) 



we obtain 



PLy{k) 



1 



[b + e'''PLy{k) + e-'''PLy{k) + 



PLy-l{k)PLy{k)]. 



(3) 



This gives 
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= — — — 7Vv 

la — Z cos k — FLy_i[k) 
which must be solved with the initial condition Po{k) = 1 for any k. For Ly < oo, the 
continued fraction expansion of Eq. (Q) must be solved numerically in general. Even in the 
zero drift case the expansion converges rapidly, as discussed in the Appendix. In the special 
case of an infinitely long terrace Ly = oo, Piyik) = pLy^i^k), and Eq. gives 



Poo{k) = a — cos k — y{a — cos k^ — b. (5) 

In Fig. 3 we show the behavior of Pl^ (x) for various values of Ly and id- In the continuum 
limit, the tail of this function goes as for the case of zero drift. 

In practice, we also need the propagator for a periodic system with a finite width L^. 
This is most easily obtained in the Fourier space by 



oo 



PLx,Ly{x)= Plyix + rL^) 

r=— oo 

= E / dke-^'^^+^'^^^PLyik) 

r=-oo ^0 

= / dke-^'^Pr.y{k)- Y: e 



oo 

ikrLx 



r=—oo 

r2i^ ~ 1 °° 97rr) 

= / dke-^'^P,y{k)±- Y 6{k-'-P) 

-'0 J^x n=-oo 

= -Y: e-'—PLy{-^)- (6) 
Numerically, Eq. (pi) is easy to implement using the Fast Fourier Transform algorithm. 



C. Continuum Limit of the Model 



It is relatively straightforward to write down a continuum description by using the diffu- 



sion equation (an electrostatic analogy |^ can also be employed). The probability density 



of a random walker u{r,t) obeys the biased diffusion equation with a source term p{f,t): 

- V ■ D ■ VM(f, t)+v- Wu{f, t) = p{f, t). (7) 
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Using the distribution of the biased random walk, we can derive expressions for the drift 
term v = (0, vj) and the diagonal elements of the diffusion tensor D = D^^, {fi, u = x,y) 



to be Vd = -id/il + id), and D,, = 1/(4 + Ma) and Dyy = 1/(4 + Ma) + 8^^/(4 + 4id 



\2 



We note that in the model, diffusion is always only slightly anisotropic for id > 0, and 
D XX /Dyy = 1/3 for £^ — > oo The source term p in Eq. (^ is constant over the whole 



terrace. The boundary conditions are that for the step edge m = 0, and for the reflecting 
boundary du/dy = 0. Also, the arrival probability of a random walker at the step edge is 
proportional to the normal derivative of the probability field u. With zero drift {id = 0), 
Eq. (0) reduces to the Poisson equation obeyed by many growth models (see e.g. Refs. 
^p6|). The present sticking rules guarantee that the growing step forms a compact 



structure, in contrast to the typical Diffusion Limited Aggregation models p3[. It is also 
evident from the stability analysis of Salditt and Spohn (see also Ref. ||ll|) that the 
one-sided diffusion field is highly destabilizing, and the step edge morphology is always 



controlled by the instability rather than described by the nonlinear KPZ equation |]T2[ . 



III. NUMERICAL RESULTS 



A. Ledge Roughness 



We have performed extensive Monte Carlo simulations of the model with the algorithm 
described in Sec. II. In this work, we consider the case of finite id only ||2^. First, we 



discuss results for the roughness of the growing ledge on an infinitely long terrace {Ly = oo) 
with a large value of Lx = 10^. In this case, after a short initial transient the undulations 
of the ledge grow and finger-like structures emerge, separated by deep cracks. The cracks 
deepen and the fingers themselves coarsen at the expense of other fingers. In Fig. 4, we 
show a sequence of typical successive configurations for different values of id- We find that 
the width w{t) of the interface associated with the ledge follows power law behavior 



w{t) = {[h{x,t) - h{t)]^y^^ = At''', (8) 
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where the brackets and the overbar denote an average over the configurations and over 
each finite system, respectively. The height variable h{x, t) is the column height of the ledge 
as measured from y = 0. Numerically, we find that the width w{t) grows linearly with 
[3i = 1.0 ± 0.01 and its slope A{id) depends on the drift id (Fig. 5). Linear growth can 
be understood qualitatively, since particles arriving at the vertical section of the ledge do 
not contribute significantly to the ledge roughness. The increase in the roughness is mainly 
due to particles that stick on top of the columns, and thus the width grows proportional to 
the total particle number, i.e. time. The change in the growth rate is also easy to explain 
qualitatively. With small drifts, only the top of the finger grows and very few particles 
reach the bottom of the cracks. With larger drifts, the probability of reaching the bottom 
increases, and thus w increases more slowly (see also Fig. 4). 

The value of /?i = 1 is consistent with the theory of Elkinani and Villain for a simple 



ID Zeno model of MBE growth with step edge barriers. Instead of ledges, they consider 
deposition of adatoms on a stepped surface with diffusion. They show that with strong step 
edge barriers, deep cracks are formed on the surface whose depth grows linearly in time. 
In this case, the deposition noise is not relevant and this result can be obtained from a 
deterministic model. 

To study the effect of a finite terrace length Ly < oo, we have simulated the model 
with id = 1/4, = 10^ and Ly = 50, 70, 100, 140, 200, and 500. Due to the fact that 
in such finite systems the relative proportion of the fiux deposited in between the fingers 
increases with time, the width w(t) eventually saturates to an Ly dependent value, but does 
not saturate as a function of L^. We find that the width satisfies the scaling ansatz of 



Family- Vicsek p,29 



w 



{Ly,t)=t^^/^^MLy/t'/^^), (9) 



where the scaling function /i(x) behaves as 



f const. X ^ 1; 
fiix) ~ (10) 
I x^'^ X <^ 1. 
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The exponent xi characterizes the surface morphology in the saturated regime w{Ly) ~ 
L^i, and the crossover time tgat ~ determines where the saturation takes over. The 
growth exponent for t <^ tsat is Pi = Xi/^i- We find that setting /5i = 1, zi = 1.00 ± 0.03 
collapses our data best to a single scaling function shown in Fig. 6. We have also obtained 
the exponent Xi by fitting to the saturated width w{Ly) and find that Xi = 0-96 ± 0.02. 

B. Finger Coarsening 

In our model diffusion along the ledge is limited by infinitely strong barriers, since the 
particles can never go around corners. This is basically the same effect as step barriers along 
the surface of the ID Zeno model. However, since in our model there is a real diffusion field 
surrounding the fingers on the terrace, additional finger coarsening takes place as is 
evident in the configurations of Fig. 4. For a finite system with < oo, this eventually 
leads to a configuration where there is only one finger present. To investigate the temporal 
scaling of the thickness of the fingers, we have studied how the first zero of the Green's 
function at r = ro{t) 

G{r, t) = {^J2 + ^> t)h{x, t) - h{tf) (11) 

X 

behaves as a function of time. The behavior of ro{t) should indicate the existence of a char- 
acteristic, time- dependent correlation length in the direction perpendicular to the direction 
of growth. In Fig. 7 we show ro(t) for several values of id when = 10^ and Ly = oo. To a 
good degree of accuracy, we find that ro(t) ~ t^'', with the value Pr = 0.32 ± 0.01 for drifts 
varying from 1/8 to 32. There is, however, a long crossover regime at the beginning of the 
growth that depends on the drift, being longer for larger drift values. 

It is also interesting to study the scaling of the Green's function. Asymptotically, we 
expect G{r, t) to scale as [^] 

G{r,t)=t-'^^ge,{t-^^r), (12) 
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where gi^{x) is a new scaling function associated with the coarsening process. In Fig. 
8 we show scahng of the data for G{r,t), with a very good data collapse obtained with 
f3i = 1 and jSr = 1/3 [0. It is interesting to note that the finger coarsening in the present 
model follows the same power law of 1/3 as Model B, which describes domain coarsening 
due to long range diffusion |^T],^. However, although qualitatively similar, the present 



scaling function depicted in Fig. 8 is quantitatively different from that of Model B |^3[. The 
exponent 1/3 also appears in models of noise-driven coarsening of mounds in ID surface 
growth where slope selection occurs because of step edge barriers [[7|p6|. 

The finite-size scaling of tq is different from that of the width w, since for a system with 
a finite terrace width < oo but with Ly = oo, the late-time configuration consist of one 
finger only, whose vertical roughness w keeps on growing linearly but whose Tq saturates. 
This introduces a new exponent Z2 that controls the saturation of tq in the x direction. On 
the other hand, for = oo and Ly < oo, both w and ro saturate, and their saturation must 
be characterized by the same exponent zi in Eq. (^. Thus, for the general case of both 
Lx,Ly < oo, we expect the following scaling form to hold: 

ro{L^,Ly,t)=tf^^M^,^). (13) 

We will not study the whole scaling function fr{x,y) here but consider the effects of a 
finite Lx and Ly separately |Q. For Ly = oo, we can again write down the Family- Viscek 



form as 

r,{L,,t)=t^^/^^h{L,/e'^^), (14) 

where now f3r = X2/Z2, and the scaling function /2 has the same limits as /i, but now 
with a new roughness exponent X2- We have simulated the model with L^ = 20, 50, 100, 
and 200 using the drift id = 1- Because of the single finger final configuration, there are 
large fiuctuations in the data and thus we have determined the saturation exponent X2 by 
estimating the saturated width rQ{Lx) directly for various values of L^. From the data, our 
best estimate is X2 = 1-02 ± 0.01, i.e. the width of the final finger grows as the horizontal 
system size. Together with f3r = 0.33 this implies that Z2 = 3.0. 

11 



In the case of a finite Ly, we expect that the scaling form satisfies 

roiLy,t)=t^^/^^h{Ly/f/^^), (15) 

where (5r must now satisfy the relation jSr = Xs/^i? with X3 being another new roughness 
exponent. The limits of and /i are again of the same form. By using system sizes Ly = 
50, 70, 100, 140, 200, and 500 with the drift id = 1/4, our data collapses to the scaling 
form shown in Fig. 9 with xs = 0.33 ± 0.01 and Zi = 1.00 ± 0.03. Moreover, we have 
obtained another estimate of the new saturation exponent xs by estimating the saturated 
width ro{Ly) and indeed verify that xs = 0.34 ± 0.02. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have in this work introduced and examined a very simple model for the 
growth of an isolated step with infinitely strong step edge barriers. The destabilizing effect 
of the one-sided biased diffusion field coupled with strongly anisotropic adatom dynamics 
makes the ledge morphologically unstable, with finger-like structures developing separated 
by deep cracks. After an initial early-time transient the fingers coarsen as t'^'^^ and the width 
of the ledge grows linearly. For an infinitely wide and long terrace, the fingers eventually 
become needle-like. We have also studied the finite-size scaling of both the coarsening and 
the width of the ledge in detail, and determined the corresponding scaling exponents. 

Recently, Pierre-Louis et al. [|14| have considered in detail a more realistic model of step 
train growth in the case of weak desorption, and one-side attachment. As in the present 
case, they find that the step morphology is linearly unstable, but now the individual step 
widths grow oc t^''^, with the steps "locked in" together. In this regime, there is no step 
coarsening, either. Thus, we expect our model to be relevant only for the case where the 
steps are well isolated, and detachment from step edges can be neglected. 

Experimentally, growth of steps on slightly miscut Si(OOl) surfaces has been studied, 
with the claimed result that the step roughening is consistent with the KPZ prediction . 
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However, at least superficially the steps depicted in Ref. [§5| appear to develop finger-like 
structures separated by deep grooves characteristic of the unstable regime studied here and 
in Ref. It would be interesting to carry out more systematic studies of roughening of 
widely spaced steps on semicondutor surfaces to characterize the nature of the instability. 
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APPENDIX: 

In this Appendix, we calculate explicitly the average arrival time tarr of a walker to 
demonstrate the need to use the present algorithm. We will also discuss the convergence of 
the probability distribution for finite terraces with Ly < oo. To begin with, the distribution 
for the number of steps or the distribution of the first passage time Piyit), can be calculated 
similarly to the spatial distribution of Eq. ([1|) by the recursion 

oo 

^'^ s=l 

oo 

J2PLy{s-u)PLy-l{u)]. (Al) 

u=l 

Using the temporal Fourier transform 

oo 

Ki^) = E e-'^'PLyit) (A2) 

t=— oo 

we obtain the characteristic function as 

pLviuj) = ^—^ , (A3) 

which can be solved with the initial value Pq{uj) = 1. Again, Piyiuj) = Piy^iiuo) when 
Ly —>■ oo and Poo(^) can be obtained. The average arrival time is proportional to the first 
derivative of the characteristic function at = by 
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u>=0\ 



b du 

when id > 0. In the infinite terrace hmit {Ly oo), iarr = 1 + while for id — it 
is easy to show that tarr = 4L2^. Thus, the return time quickly becomes prohibitively large 
for large systems and small values of the drift, making brute- force Monte Carlo simulations 
difficult. On the other hand, for drifts larger than unity, no significant reduction in computer 
time can be obtained with the new algorithm. 

Finally, to estimate the convergence of the probability distribution PLy{x) towards its 
asymptotic limit as a function of the terrace length Ly for any L^, we can define the deviation 
d by 

d {Lx, Ly, id) 

- r 't\PLyil^) - ^oo(^)]^ (A5) 

assuming a periodic system in the x direction. In Fig. 10, we show the deviation for 
various values of the drift as a function of Ly. 



14 



REFERENCES 



[1] D. J. Eaglesham, J. Appl. Phys. 77, 3597 (1995). 

[2] A.-L. Barabasi and H. E. Stanley, Fractal Concepts in Surface Growth, (Cambridge 
University Press, 1995). 

[3] Reflection High Energy Electron Diffraction and Reflection Electron Imaging of Surfaces, 
ed. by P. K. Larsen and P. J. Dobson (Plenum Press, New York, 1988) 

[4] W. K. Burton, N. Cabrera, and F. C. Prank, Phil. Trans. R. Soc. London Ser. A 243, 
299 (1951). 

[5] R. Chez and S. S. Iyer, IBM J. Res. Develop. 32, 804 (1988). 
[6] R. L. Schwoebel, J. Appl. Phys. 44, 614 (1969). 

[7] J. Krug, M. Phschke, and M. Siegert, Phys. Rev. Lett. 70, 3271 (1993); J. Krug and H. 
Dobbs, ibid. 73, 1947 (1994); M. Siegert and M. Plischke, Phys. Rev. E 50, 917 (1994). 

[8] M. Siegert and M. Phschke, Phys. Rev. Lett. 73, 1517 (1994). 

[9] M. Kotrla and P. Smilauer, Acta Phys. Slovaca 44, 237 (1994); F. Liu and H. Metiu, 
Phys. Rev. B 49, 2601 (1994); C. J. Lanczycki and S. Das Sarma, Phys. Rev. E 50, 213 
(1994); S. A. Chalmers, J. Y. Tsao, and A. C. Gossard, J. Appl. Phys. 73, 7351 (1993); 
D. E. Wolf, Phys. Rev. Lett. 67, 1783 (1991); L.-H. Tang and T. Nattermann, Phys. 
Rev. Lett. 66, 2899 (1991); A. K. Myers-Beaghton and D. D. Vvedensky, Phys. Rev. B 
42, 9720 (1990). 

[10] I. Bena, C. Misbah, and A. Valence, Phys. Rev. B 47, 7408 (1993); A. Karma and C. 
Misbah, Phys. Rev. Lett. 71, 3810 (1993). 

[11] G. S. Bales and A. Zangwill, Phys. Rev. B 41, 5500 (1990). 

[12] T. Salditt and H. Spohn, Phys. Rev. E 47, 3524 (1993). 



15 



[13] O. Pierre- Louis and C. Misbah, Phys. Rev. Lett. 76, 4761 (1996). 

[14] O. Pierre- Louis, C. Misbah, Y. Saito, J. Krug, and P. Politi, unpublished (1997). 

[15] M. Kardar, G. Parisi, and Y.-C. Zhang, Phys. Rev. Lett. 56, 889 (1986). 

[16] H. J. W. Zandvliet and H. B. Elswijk, Phys. Rev. B 48, 14269 (1993). 

[17] In the case of one-particle diffusion, the drift term is not strictly speaking present since 
the step does not move. In our model £d may be considered as a finite-density correction 
to this limit. 

[18] C. Roland and G. H. Gilmer, Phys. Rev. Lett. 67, 3188 (1991); Phys. Rev. B 46, 13428 
(1992); D. Srivastava and B. Garrison, Phys. Rev. B 47, 4464 (1993); C. P. Toh, S. H. 
Seow, and C. K. Ong, Surf. Sci. 292, 114 (1993). 

[19] Q.-M. Zhang, C. Roland, P. Boguslawski, and J. Bernhole, Phys. Rev. Lett. 75, 101 
(1995). 

[20] J. Hoeven, J. M. Lenssinck, D. Dijkkamp, E. J. van Loenen, and J. Dieleman, Phys. 
Rev. Lett. 63, 1830 (1989); Y.-W. Mo, B. S. Swartzentruber, R. Kariotis, M. B. Webb, 
and M. G. Lagally, Phys. Rev. Lett. 63, 2393 (1989); Y. W. Mo, J. Kleiner, M. B. Webb, 
and M. G. Lagally, Phys. Rev. Lett. 66, 1998 (1991); Y.-W. Mo and M. G. Lagally, 
Surf. Sci. 248, 313 (1991); Y.-W. Mo, J. Kleiner, M. B. Webb, and M. G. Lagally, Surf. 
Sci. 268, 275 (1992); B. Voigtlander, T. Weber, P. Smilauer, and D. E. Wolf, Phys. 
Rev. Lett. 78, 2164 (1997). 

[21] J. E. Griffith, G. P. Kochanski, J. A. Gubby, and P. E. Wierenga, J. Vac. Sci. Technol. 
A 7, 1914 (1989). 

[22] G. Brocks, P. J. Kelly, and R. Car, Phys. Rev. Lett. 66, 1729 (1991); C. P. Toh and C. 
K. Ong, Phys. Rev. B 45, 11120 (1992). 

[23] T. A. Witten and L. M. Sander, Phys. Rev. Lett. 47, 1400 (1981). 



16 



[24] On a Si(OOl) terrace diffusion along the dimer rows is much faster than across them. 
For ifi > 0, however, the step morphology is controlled by the drift term, and thus the 
exact value of the ratio D^x/Dyy is not important. 

[25] J. Krug, K. Kassner, P. Meakin, and F. Family, Europhys. Lett. 24, 527 (1993). 

J. Krug, Adv. Phys. 46, 139 (1997). 



[27] The zero drift limit — is rather subtle in the present model. For systems with finite 
Lx, we find that asymptotically w{t) ~ t. However, the behavior of the finger coarsening 
is somewhat different from the case of id > 0. We do find an intermediate time regime 
where ro{t) ~ but for larger systems there seems to be a late-time crossover to a 
larger effective coarsening exponent. In this regime, our model becomes equivalent to 
the "Laplacian needle model A" of Krug et al. |^5|]. In the needle model, the average 



distance between the needles grows like ha/ \og{ha), where ha is the height of "active" 



needles [Q. This would mean that ro(t) ~ t/log(t), which is qualitatively supported 
by our simulation results for i^ = 0. 

[28] I. Elkinani and J. Villain, J. Phys. I France 4, 949 (1994). 

[29] F. Family and T. Vicsek, J. Phys. A 18, L75 (1985). 

[30] P. Politi and J. Villain, Phys. Rev. B 54, 5114 (1996). 

[31] A. J. Bray, Adv. Phys. 43, 357 (1994). 

[32] As defined in Eq. (13), the scaling function ge^{x) depends on the drift id- We can also 
write the scaling form as G{r, t) = w^{id, t)g{r /rQ^idi iii which case the scaling plots 
for the new function g{x) collapse within the accuracy of our data for the values of the 
drift studied, namely Id = 1/8, 1/4, 1/2, and 1. 

[33] T. M. Rogers, K. R. Elder, and R. C. Desai, Phys. Rev. B 37, 9638 (1988). 

[34] In the long-time limit, for an infinite system [L^ = Ly = oo) the aspect ratio 



17 



tq/w approaches zero, and the fingers become needle-hke even for finite values of 
^d■ In the needle model, the penetration probability of the diffusing particles be- 
tween needles decays exponentially as exp(— vrr/^) where the distance r is mea- 
sured from the top of the needles and ^ is the distance between two needles. In our 
model with a finite drift parameter ^^j, the corresponding result can be shown to be 



exp{-2^[W^^|2^^ + where Vd{t), D^^, and Dyy are defined in 



the text. For zero drift (where Vd = 0, and Dyy = D^x), this correctly reduces to the 
previous result. 



[35] F. Wu, S. G. Jaloviar, D. E. Savage, and M. G. Lagally, Phys. Rev. Lett. 71, 4190 
(1993). 




18 



FIGURES 

FIG. 1. A schematic view on adatom dynamics on a vicinal surface. The local adatom 
concentration is denoted by n and the other symbols are explained in text. 

FIG. 2. Adatom dynamics in the growth model. For Ly < oo, deposition occurs uniformly 
randomly at all available (unoccupied) sites, while for Ly = oo, the particles are released from the 
line at y = H + 1. Particles then diffuse on the lower terrace, and drift in the —y direction. They 
become incorporated into the step when they either land on a top section of the step (along the x 
axis), or slide down along the —y direction to the nearest kink site, as shown schematically in the 
figure. If a particle attempts to cross the line at {xo,H + 1) in the y direction, it is immediately 
returned from (xq, .ff + 2) to the line with new coordinates (xq + x, if + 1) chosen from the spatial 
distribution PLy{x). The boundary above y = H + Ly + 1 (not shown) is completely reflecting. 

FIG. 3. Probability distribution Poo{x) for the return position with L^ = 10^ and Ly = oo. 
The drift parameter id = 1, 1/2, 1/4, 1/8, 1/16, and (from top to bottom at a; = 0). In the limit 
— > oo the distribution approaches a delta function. 

FIG. 4. Ten consecutive step profiles from the growth model with L^ = 10^ and = oo at 
t = 1000,2000, 10"^ with the drift (a) id = 1/4, (b) id = 1, and (c) id = 4. Only part of the 
system is shown. 

FIG. 5. The step width w{t) for 4 = 1/8, 1/4, 1/2, 1, 2, 4, 8, 16, and 32 (from top to bottom) 
with Lx = 10^ and Ly = oo. The slope A as a function of the drift id is shown in the inset. 

FIG. 6. Scaling function /i of Eq. (10) for the step width w{Ly, t) with Ly = 50, 70, 100, 140, 
200, and 500. Good scaling is obtained with Pi = 1 and z\ = l. The drift = 1/4 and the lateral 
lattice size Lx = 10^. 

FIG. 7. The finger width ro(t) for id = 1/8, 1/4, 1/2, 1, 2, 4, 8, 16, and 32 (from top to 
bottom) with Lx = 10^ and Ly = oo. The dashed line indicates a slope of 1/3. 
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FIG. 8. Scaling function g^^^ of Eq. (13) for G{r,t) with = 10^ and = ^ ten different 
times t = 1000, 2000, ...,10^, with /?i = 1 and l3r = 1/3. 

FIG. 9. Scaling function /a of Eq. (16) for the finger width vq with the terrace length Ly = 50, 
70, 100, 140, 200, and 500. The drift 4 = 1/4 and the lateral lattice size = 10^, and = 1/3 
and zi = 1. 

FIG. 10. Root mean square deviation d of the propagator Pi^y from P^o with = 10^ for 
4 = 0, 1/16, 1/8, 1/4, and 4 (from top to bottom), shown as a function of the distance Ly to the 
reflecting boundary. 
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